Realization of decoherence-free subspace using Multiple-Quantum coherences 
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This letter presents a two-dimensional nuclear magnetic resonance(NMR) approach for construct- 
ing a two-logical-qubit decoherence-free subspace (DFS) based on the fact that the three protons 
in a CH3 spin system can not be resolved in one-dimension NMR spectroscopy, but to a certain 
extent, can be distinguished by two-dimensional multiple-quantum NMR. We used four noisy phys- 
ical nuclear spins, including three protons and one carbon in the CH3 spin system, to generate two 
decoherence-free logical qubits. It made full use of the unaddressed spins which could not be used in 
one-dimensional spectrum. Furthermore, we have experimentally demonstrated such an approach. 
Our experimental results have shown that our DFS can protect against far more types of decoher- 
ence than the one composed of four noisy physical qubits all with different chemical shifts. More 
importantly, this idea may provide new insights into extending qubit systems in the sense that it 
effectively utilizes the magnetically equivalent nuclei. 
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Quantum computers could, in theory, be superior to classical computers when performing a number of computational 
tasks, such as factorizing a large number |l|, searching unsorted database |2j and, especially, simulating quantum 
systems themselves Q- Among various techniques, liquid nuclear magnetic resonance (NMR) approach is, currently, 
the most realistic tool Q| . If the chemical shift of a nuclear spin is well resolved and can be controlled by selective radio- 
frequency (RF) pulses without disturbing the other spins, the spin can be used as a qubit for quantum computing. 
Although NMR quantum computing has made the greatest progress towards QIP [4J, it is very difficult to scale up 
current NMR quantum computers to large sizes (more than 10 qubits) 0]. There are many obstacles that limit wide 
applications of NMR in quantum computation. Firstly, in addition to the low signal-to-noise ratio (SNR) in NMR 
experiments, the SNR decreases exponentially with increasing of qubit numbers. Secondly, there are only few types of 
nuclear spins which might be used as qubits, such as 1 H, 13 C, 15 N, 19 F, and 31 P etc. Thirdly, because of limitations in 
the number of the receive channels and the transmit channels in the commercial NMR spectrometers, one always had 
to employ homonuclear systems and selectively control each of the spins when implementing multi-qubit quantum 
computing. Finally, frequency dispersion limits the number of usable homonuclear spins or qubits. In the case of 
proton NMR at 500 MHz, the chemical shift range is generally less than 5000 Hz. This small spectral window restricts 
the number of qubits in a homonuclear system. In addition, it is extremely hard to select or design a special molecule 
whose protons all have well resolved chemical shifts and are coupled to each other. It is thus a great challenge for 
NMR spectroscopists to utilize all of the available spins in the real molecule. It has been noticed that magnetically 
equivalent spins, such as the three protons in a methyl, are used as a single physical qubit. One can wonder whether 
it is possible to use the equivalent spins as individual qubits. In this context, we will describe a two-dimensional (2D) 
NMR technique to implement logical qubit by making full use of the three protons in the methyl. 

The new concept of using 2D NMR to perform quantum computing was first proposed by Madi and Ernst et al. 
in 1998 6J. The quantum logical gates and Deutsch-Jozsa quantum algorithm were also experimentally realized in 
2D NMR with the help of spin- and transition- selective pulses0,Q- In these experiments, since single quantum 
transitions were used in both dimensions, the equivalent protons are indistinguishable and can only be utilized as 
one physical qubit. On the contrary, using 2D multiple-quantum NMR J-resolved spectroscopy (MQ-JRES) of the 
spins in alanine 9], we observed that the methyl has four resolved multi-quantum (MQ) coherences which have no 
interchange during the free evolution time. The 2D MQ-JRES spectrum can be explained in detail by a SI3-M spin 
system (for the 13 CH3- 12 CH in alanine). After four spins are excited to their highest quantum state.the SI3 spin 
system will include four coherences (one quadruple quantum, two double quanta, and one zero quantum) 9]. Moreover, 
different multi-quantum coherences are modulated differently by their coupling constants to the remote spin M, and 
is distinguished by projecting the 2D multi-quantum onto the Fl dimension. In this case, the three protons in the 
methyl play the roles of three physical qubits. The following question arises: "can three physical qubits composed of 
three protons in a methyl be used to implement the 2D NMR quantum computing(QC)?" This topic will be discussed 
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and solved in the following context. 

Here, we will describe how to use three protons and a carbon in a methyl as four physical qubits to produce four 
multiple quantum coherences, and then utilize these multiple quantum coherences as logical qubits to construct an 
operator subspace DFS in a Liouville space0, ||| ^| by virtue of the 2D multi-quantum J-resolved NMR spectrum 
(MQ-JRES). We also present experimental results in demonstrating this novel idea. 

There are two kinds of passive error control codes for counteracting the decoherence. One is decoherence-free 
subspace (DFS) 0,0,0 and the other is noiseless subsystems (NSs) [l4 H5I Hil Il7| . They had been widely studied 
in both theoretically and experimentally [ill, |l2|, lljj, LLJ, Ilia llo, LLa ■ Recently the two-logic-qubit DFS composed 
of four physical qubits were demonstrated in NMR quantum computing^]. All states in DFS are immune to the 
error operation E,j with the form Ed — otdfiExE^E^Ei + ad, 1X1X2 E3E4 + ad^EiE^X^Xi + 04,3X1X2X3X4, where E 
is the unit operator and Xi (i=l, 2, 3, and 4) are operators which make the physical qubits flipped. Just as the other 
methods to avoid the decoherence problem, the construction of the DFS always needs at least two physical qubits to 
produce one logical qubit. 

In order to express a DFS as state vectors or density matrices, it is instructive to discuss the connection between 
logical-qubit state vectors and density matrices. For simplicity, we first consider a DFS with two physical qubits jj^, 
which protects against the error set {E1E2, X1X2, Z1Z2} in the form of 



|o>£ = _ (|oi) + |io), 

|1>£ = ^(|01)-|10). 



(1) 



The above states can also be rewritten in the form of density operators and the state |0)° corresponds to the following 
expression, 



fa) = lo>M 



(2) 



On the other hand, an experimental density matrix (p) could be separated into different orders p ( such as zero-, 
single-, two-, or quadruple-quantum coherence and so on)|2lj. 

+21 

p= E * 

p=-2l 

where I is the maximum total spin quantum number. For example, we can decompose Eq. 2 into the superposition 
of zero quantum and double quantum, 



n L - n° -I- n 2 -I- - - T 1 T 2 +- T 1 T 2 -I- T 1 T 2 

P(L0) — P +P +P — J ~ 1 z 1 z + 1 x 1 x + 1 y 1 y 

P 2 +p- 



(3) 



where I z , I x , I y are product operators with I z =-^a z 
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In a 



NMR experiment the different orders will evolve according to their own rules. Actually, the different orders are 
the elementary constituents of the density matrices. In this case, it is possible to create a DFS just by using the 
different orders as basis states in the density operator representation. Therefore, we should first find a group of logical 
orthonormal basis states which span a self-contained space. Such a group of basis states may be expressed in the form 
of density operators just as the four logical states proposed in Ref.[18|. But, on the other hand, it could also exist in 
an operator subspace and be presented in the fashion of density operators as long as one can find a group of logical 
orthonormal basis states. 

In looking for such an operator subspace DFS, we noted that there existed four MQ coherences in the spin system 
SI3 (such as the methyl 13 CH3 where S denotes 13 C and I stands for 1 H), including one quadruple quantum (QQ), 
two double quantum (DQ1,DQ2) and one zero quantum (ZQ). More importantly, these four MQ coherences can be 
distinguished by applying Liu et aVs pulse sequence^]. Their product operator forms are 
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Pl(QQ) — 3I x IyIySy 3I x IxIySx IxIxIxSy ~f~ lylylySxi 

P2(DQ1) = 3I x IyIySy ~\- ?)I x IxIySx I X Ix^xSy IyIyIyS x , 

P3(DQ2) — Ix^yly^y ^-x^-x^ySx ~t~ ^x^x-^x^y ~l~ lylylySx: 
P4(^Q) ^x-^yly^y Ixlxly^x ~t~ Ix^-x^x^y -^y^y^y^x- (4) 

It is obvious that not all of the above four MQ coherences are orthogonal to each other, namely, they could not be 
used as orthonormal basis states for a four-dimension operator space. Nevertheless, if we slightly modify them to the 
following set of operators shown in Eq. 5, it is interesting to find that the four modified density operators (they also 
belong to multiple quantum coherences) span a two-logical-qubit DFS. We denote the four modified MQ coherences 
as p 1} p 2 , p 3 , P4, which correspond to the logical states |00)ll(00|, |01)ll(01|, |10)ll(10|, respectively. 
These states are written as 



Pl = |00)LL(00| = I X IyIySy — I X IxIyS X — IxIxIxSy + IyIyI y S X , 
P2 = |01)LL(01| = I X IyIySy + IxIxIySx ~ ^X^X^xSy — IyIyIyS X , 

Pa = |10)ll(10| — Ix ^y ^y Sy ~\~ I x Ix ly S x -\- I x Ix Ix Sy ly ly ly S x ) 

P4 = |H)ll(H| = IxIylySy I X I X IyS X + I X IxIxSy IyIyIyS X , (5) 

where the normalization constants have been omitted. Specifically, the states corresponding to these density operators 
are the eigenstates of the operator E n , 



E n = a nfi E 1 E 2 E 3 E i +a nA E 1 E 2 E 3 Z 4 +a, h2 Z 1 Z 2 Z 3 E i +a n - i Z 1 Z 2 Z 3 Z i + 

a nA X 1 X 2 X 3 X 4 +a n . 5 X 1 X 2 X 3 Y i +a nfi Y 1 Y 2 Y 3 X 4 +a n jY 1 Y 2 Y 3 Y i . (6) 

In a word, by using the spin system SI3 we have actually constructed a two-logical-qubit DFS which can resist the 
error E n . Now the next work is to find an effective experimental technique of realizing the DFS with the above 
product operators. 

Our method to implement the above DFS is based on the idea in Ref . || , in which the pulse sequence can be divided 
into three parts. The first part (the pulse sequences before the encoding gradient pulse G e ) is to generate the highest 
state 8I x IyI y S y which consists of the above four coherences. In the second part (including the gradient pulse G e , the 
evolution time ti and the two 180° pulses on spins / and S), each of the above four coherences evolves differently. At 
the end of the second part, the four coherences have the following forms 



Pl(QQ) {^IxIylySy 3I X I X IyS X Ix^IxSy ~\~ ly ly ly S X ^j 

cos[7r(3J/ju + J S M)h] exp(-ti/T 2QQ ), 

P2(-DQ1) = {3I X IyIySy + 3I X I X IyS X IxIxIxSy IyIyIyS X ) 

cos[tt(3J/m - Jsu)h] exp(-ti/T 2 uQi), 

P3(DQ2) — {IxIylySy + I x IxIyS x + I x I x I x S y + I y I y I y S x ) 

cos[7r( Jim + JsAi)h} exp(-ti/T 2DQ2 ), 

P4:(ZQ) — {IxIylySy I x lxly$x H - Ix^-x^-xSy lylylySx} 

cos[7r( J/m - JsmYx] exp(-£i/T 2 zQ). (7) 

Therefore, the four MQ coherences can be distinguished by the remote J couplings of the spins I and S with the spin 
M, respectively, and the transverse relaxations. The third part is to transform the multiple-quantum coherences into 
the detectable signals. It should be noted that different multiple-quantum coherences have different projection onto 
the Fi dimension, which means that the resonances of the different multiple-quantum are resolved. In particular, 
the selectively detected MQ-JRES spectrum can be achieved by adjusting the strength ratio between the encoding 
and decoding gradient pulses(i.e. between G e and G c i)- One example is that when the gradient strength ratio ^ is 
-8:10, the double quantum coherence(DQ2) could be selectively detected. It is also worthwhile noting that there is 
no interchange between the multi-quantum coherences during the evolution period. 
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Considering that product operators p i in Eq. (4) and pj(MQ) in Eq. (5) have the similar forms, the states 
p 1 ,p 2 , / o 3 ,P4 could be obtained by using special pulse sequences similar to the one in Rcf. 9]. That is, the method 
used to yield the states p 3 and p 4 is, respectively, the same as the one used to produce the multi-quantum coherences 
p 3 (DQ2) and p 4 (ZQ). The states p x and p 2 can also be created by applying the pulse sequences corresponding to 
the multi-quantum coherences p x (QQ) and p 2 (DQl) apart from adding two special-angle pulses on spins C and H, 
respectively. 

Interestingly, due to the same chemical shift of the three protons in a methyl, error operators only affect them 
simultaneously, and the following error operator E m doesn't exist: 

E m = a m , (T 1 E 2 E 3 E 4 +a mA E 1 c7 2 E 3 E 4 +a m ^E 1 E 2 a 3 E 4 +a m ^a 1 a 2 E 3 E 4 + 
a mA a 1 E 2 cj a E 4 +a m ^E 1 a 2 c7 3 E 4 +a mfi (J 1 E 2 E 3 a 4 +a m ^E 1 E 2 c7 3 a 4 

, pi 2p3 4 1 S t?3 A , 1 t-,2 3 A , j^l 2 3 A /o\ 

+a m ^E a E a -\-a m $cr a E a +a m ^a E a a +a m> iiE a a a , (8) 

where a, a, a S {X, Y, Z}. In other words, the DFS we constructed can avoid errors E n and E m simultaneously. 
However, that one proposed by Ollerenshaw et al. based on four different physical qubits only protests against small 
parts of errors E n and E m . It implies that it is more effective to adopt the SI n system to realize the DFS. On the 
other hand, in ID NMR quantum computing, the spins I n in SI n (n=2, 3) were always regarded as one qubit, so in a 
sense we actually implemented a DFS in a two-qubit physical quantum system, which has more advantages over the 
one constructed by four physical qubits. 

The experiment was carried out on a Varian INOVA 600 spectrometer. The sample was 20mg alanine in 0.5ml D2O, 
with one carbon in methyl labelled. The molecular structure is 13 CH3- 12 CH(NH 2 )- 12 COOII. The spin system SI 3 M 
presents 13 CH3- 12 CH in the alanine molecule. The coupling constants obtained from the high-resolution ID NMR 
spectra, Jsi, Jsm and Jim are 129.8Hz, 4.5Hz and 7.3Hz, respectively. Using pulsed field gradients, which is similar to 
the method in Ref.|I|, we can selectively obtain the states p 1 , p 2 , p 3 or p 4 . For example, p 3 can be selectively detected 
by applying the pulse sequences shown in Fig.l. Here, the required gradient strength ratio between the encoding 
gradient(G e ) and the decoding gradient (d) is -8:10. The delay A is 3.88ms corresponding to the coupling constants 
Jsi=129.8Hz. The spectral widths are 30 and 4000 Hz in the Fi and F 2 dimensions, respectively. The data points in 
the Fi and F 2 dimensions are 64 and 10000. In order to show that this DFS' protection against the errors E n and 
E m , we performed a series of errors belonging to E n and E m at the position c in Fig.l. Whether the DFS constructed 
by two-dimensional MQ J-resolved spectrum protects against errors E n and E m or not, is demonstrated by observing 
and comparing the 2D NMR spectra before and after applying the decoherence operators. 

Typical experimental results are shown in Fig. 2, where (A), (B) correspond to the spectra of p 3 before and after 
one performs the error operator X 1 X 2 X 3 Y 4 , respectively. It is interesting to note that the spectrum in (A) is the 
same as the one in (B). This means that p 3 is certainly not affected by this decoherence. Similar outcomes are also 
obtained when the other operators belonging to E n or E m are applied to p l , p 2 , p 3 or p 4 . Due to limitation in length, 
we do not reproduce the experimental results. 

Decoherence has been one of the severe impediments in quantum information processing. Many techniques have 
been widely studied. Therefore it is significative to experimentally and theoretically focus on the DFS that has been 
developed as a tool for protecting QIP. The idea for constructing the DFS provided in this paper shows new insights 
into effectively utilizing the spin system including CH n (n=2,3) to achieve quantum computing as well as extending 
the number of qubits. This process is impossible in ID NMR, since the magnetically equivalent spins in CH n can 
not been fully used. But here we showed that 2D MQ-JRES could do it by using the MQ coherences which were 
modulated by the remote spin. In addition, a novel point is that since the n numbers of H have the same chemical 
shift, the errors only affect them together, so to this extent, this DFS may avoid more errors than others that were 
constructed by spins with different chemical shifts. 

Simultaneously, when one performs ID NMR quantum computing with a sample containing the CH n group, all the 
protons in each CH n group which could only be used as one qubit, couple to the other spins. In this condition, the 
spectra of the other spins become more complex than those for which there are no magnetically equivalent spins in the 
sample, although their number of qubits is the same. So it seems that the CH n group just makes ID NMR QC more 
complex. However, just as described above, the CH n group, in 2D MQJES, makes our DFS more powerful. Here, we 
merely performed a special DFS applying 2D MQJES. The demonstration of some quantum algorithms based on this 
DFS is under way. 

In summary, we have provided a new two-logical-qubit DFS with 2D MQ-JRES and the three magnetically equiv- 
alent spins in CH3. In a sense, the three protons are distinguishable. The corresponding experimental demonstration 
has been completed. Compared to the other schemes for implementing DFS, our method protected against more 
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error operators just by employing a smaller nuclear system. Specifically, this idea can also be applied to implement 
other QIP tasks on NMR QC. It follows that when the selected spin system contains a methyl, we can implement 
an (n+2)-qubit quantum computing just using n qubits of one-dimensional NMR QC. Conclusively, the 2D NMR 
method presented in this letter may be helpful for developing the multi-qubit NMR QC. 

This work has been supported by the National Natural Science Foundation of China under Grant No. 10374103, 
the National Basic Research Programme of China under Grant No.2001CB309309, and the National Knowledge 
Innovation Program of the Chinese Academy of Sciences under Grant No: KJCX2-W1-3. 
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FIG. 1: Pulse sequence for multi-quantum J-resolved NMR spectroscopy showing p 3 can protect against the error X 1 X 2 X 3 Y 4 , 
where a is the encoding process, b is the decoding process, and c is the position to perform errors. 
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FIG. 2: Multi-quantum NMR spectra of p 3 before (A) and after (B) performing the error X*X 2 X 3 Y 4 . The projection into F2 
dimension presents the 1 H NMR spectrum of the 13 CH3 group protons with the signal from 12 Cfi3 protons. The projection 
into Fi dimension shows the MQ spectrum, which is given at the right side of the 2D plot. 
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